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Abstract 

We study the asymptotic behavior of linear evolution equations of the 
type dtg = Dg + Cg — \g, where £ is the fragmentation operator, D 
is a differential operator, and A is the largest eigenvalue of the operator 
Dg + Cg. In the case Dg = —d x g, this equation is a rescaling of the 
growth-fragmentation equation, a model for cellular growth; in the case 
Dg = —d x (xg), it is known that A = 1 and the equation is the self- 
similar fragmentation equation, closely related to the self-similar behavior 
of solutions of the fragmentation equation dtf = Cf. 

By means of entropy-entropy dissipation inequalities, we give general 
conditions for g to converge exponentially fast to the steady state G of the 
linear evolution equation, suitably normalized. In other words, the linear 
operator has a spectral gap in the natural L 2 space associated to the 
steady state. We extend this spectral gap to larger spaces using a recent 
technique based on a decomposition of the operator in a dissipative part 
and a regularizing part. 
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1 Introduction and main results 

In this work we study equations which include a differential term and a fragmen- 
tation term. These equations are classical models in biology for the evolution 
of a population of cells, in polymer physics for the size distribution of poly- 
mers, and arise in other contexts where there is an interplay between growth 
and fragmentation phenomena. The literature on concrete applications is quite 
large and we refer the reader to [7, 12] as general sources on the topic, and to the 
references cited in [3, 8, 13] for particular applications. We deal with equations 
of the following type: 



The unknown is a function gt(x) which depends on the time t > and on x > 0, 
and for which an initial condition gi n is given at time t = 0. The quantity 
gt{x) represents the density of the objects under study (cells or polymers) of 
size i at a given time t. The function a = a(x) > is the growth rate of 
cells of size x. Later we will focus on the growth-fragmentation and the self- 
similar fragmentation equations, which correspond to a(x) = 1 and a(x) = x, 
respectively. 

Most importantly, we pick A to be the largest eigenvalue of the operator 
g i y —d x (ag) + Cg, acting on a function g = g(x) depending only on x; see 
below for known properties of this eigenvalue and its corresponding eigenvector. 



dtgt{x) + d x (a(x)g t {x)) + Xg t (x) = C[g t ]{x) 
9t(0) = (t > 0) 
9o(x) = gin(x) (x > 0). 



(la) 
(lb) 
(1c) 
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The fragmentation operator C acts on a function g = g(x) as 

Cg(x) := C+g(x) - B(x)g(x), (2) 
where the positive part is given by 

£+9{x) ■= / Hy,x)g(y)dy. (3) 

The coefficient b(y,x), defined for y > x > 0, is the fragmentation coefficient, 
and -E>(x) is the toiaZ fragmentation rate of cells of size x > 0. It is obtained 
from 6 through 

B(x) := / ^b(x 7 y)dy (x>0). (4) 

Asymptotic behavior As said above, we pick A to be the largest eigenvalue 
of the operator g H > —d x {ag) + Cg. Under general conditions on b and a, it is 
known [3, 8] that A is positive and has a unique associated eigenvector G (up 
to a factor, of course), which in addition is nonncgativc; i.e., there is a unique 
G solution of 

(a(x)G(x))' + XG(x) = C(G)(x) (5a) 
0(1)0(1)1^ = 0, (5b) 

G>0, / G(x)dx = l. (5c) 



The associated dual eigcnproblem reads 

-a{x)d x (f> + (B(x) + A) (j){x) = C* + cf>(x), (6a) 
<p>0, / G{x)(j){x) dx = 1, (6b) 



where 

C*+4>(x) := f b(x,y)0(y)dy, (7) 
Jo 

and we have chosen the normalization J Gtj) = 1. This dual eigenproblcm is 
interesting because <f> gives a conservation law for (1): 

/>OG 

<j>{x) gt(x) dx = I 4>{x) gi„(x) dx = Cst (t > 0). (8) 
Jo 

The eigenvector G is an equilibrium of equation (1) (a solution which does 
not depend on time) and one expects that the asymptotic behavior of (1) be 
described by this particular solution, in the sense that 

g t — !> G as t —> oo, (9) 
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with convergence understood in some sense to be specified, and g normalized 
so that / 4>(x)g t (x) dx = 1. Furthermore, one expects the above convergence to 
occur exponentially fast in time; this is, 

\\g t -G\\<C\\g in -G\\e-^, (10) 

for some (3 > 0, in some suitable norm. This latter result has been proved in 
some particular cases which are essentially limited to the case of B constant 
[13, 6]. In this paper we want to prove this result for more general B, which we 
do by using entropy methods. In order to describe our results, we need first to 
describe the entropy functional for this type of equations. 



Entropy The following general relative entropy principle [9, 10] applies to 
solutions of (1): 



d_ 
dt 



OO pOO 



{x)G{x)H (u{x)) dx = / (/)(y)b(x,y)G(x) 

JO Jy 



ID 



x (H(u(x)) - H(u(y)) + H'{u{x)) (u(y) - u(x)) ) dxdy, ( 
where H is any function and where here and below 

u{x) := |H (x> 0). (12) 

When H is a convex function, the right hand side of equation (11) is nonpositivc, 
and in the particular case of H(x) := (x — l) 2 we have 

j t H 2 [g\G] = -D b [g\G] < 0, (13) 

where we define 

poo 

H 2 [g\G}:= 4>G{u- l) 2 dx and (14) 
Jo 

/•OO />OG 

D b [g\G}:= / <f>(x) G(y) b(y, x) (u(x) - u{y)f dydx. (15) 



J x 



Since 



/*°° /A 

H 2 [g\G] = (g-G) 2 ^dx=\\g-G\\ 2 L2 ^ G - ldx) 



(16) 



we see that proving H2{gt\G) — > implies that the long time trend to equilibrium 
(9) holds, and proving the entropy-dissipation entropy inequality 

H 2 [g\G]<^D b [g\G], (17) 

implies that the exponentially fast long time trend to equilibrium (10) holds for 
the same /3 > and for the norm || • |[ = || • \\L 2 (<j>G- l dx)- 
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The main purpose of our work is precisely to study the functional inequality 
(17) and establish it under certain conditions on the fragmentation coefficient b. 
We notice that, while some results on convergence to equilibrium for equation 
(1) are available, no inequalities like (17) were known, so one of the main points 
of our work is to show that the entropy method is applicable, in certain cases, 
to give a rate of convergence for this type of equations. 

Wc focus on the two remarkable cases a(x) = 1, which corresponds to the so- 
called growth-fragmentation equation, and a(x) = x, which gives the self-similar 
fragmentation equation. There are several reasons for restricting our attention 
to them, the main one being that they are the ones most extensively studied in 
the literature due to their application as models in physics and biology. Also, 
the inequality (17) depends on the particular properties of the solution 4> to 
the dual eigenproblem (6) and the equilibrium G. The existence and properties 
of these have been studied mainly for the above two particular cases, and are 
questions that require different techniques and deserve a separate study. One 
of our main results gives general conditions under which the entropy-entropy 
dissipation inequality (17) holds, and this may be applied to cases with a more 
general a(x), once suitable bounds are proved for the corresponding profiles (j> 
and G. 

In order to understand our two model cases, let us describe them in more 
detail and give a short review of previously known results for them. 

The growth-fragmentation equation The growth-fragmentation equation 
is the following [6]: 

d t n t + d x n t = Cn tl (18a) 

M°) = (t > o) ( 18b ) 

n (x) = n in {x) (x > 0). (18c) 

Here, rit {x) represents the number density of cells of a certain size x > at time 
t > 0. The nonnegative function m n is the initial distribution of cells at time 
t = 0. Equation (18) models a set of cells which grow at a constant rate given by 
the drift term d x n, and which can break into any number of pieces, as modeled 
by the right hand side of the equation. The quantity b(x,y), for x > y > 0, 
represents the mean number of cells of size y obtained from the breakup of a 
cell of size x. If one looks for solutions of (18a)-(18b) which are of the special 
form n t (x) = G(x)e xt , for some A e R, one is led to the eigenvalue problem 
for the operator — d x + C and its dual eigenvalue problem, given by equations 
(5) and (6), respectively. It has been proved [3, 8] that quite generally there 
exists a solution to this eigenvalue problem which furthermore satisfies A > 
and G > 0. 

For this particular A, wc consider the change 

g t (x):=n t (x)e- xt . (19) 

Then, n satisfies the rescaled growth-fragmentation equation, which is the partic- 
ular case of eq. (1) with a{x) = 1. The long time convergence (9) or (10) means 
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here that the generic solutions asymptotically behave like the first eigenfunction 
G(x) e xt (with same <f> moment). 

Let us give a short review of existing results on the asymptotic behavior 
of equation (1) for constant a(x). In [9, 10] the general entropy structure of 
equation (1) and related models was studied, and used to show a result of 
convergence to the equilibrium G without rate for quite general coefficients 
b, but under the condition that the initial condition be bounded by a constant 
multiple of G [10, Theorem 4.3]. Some results were obtained later for the case of 
mitosis with a constant total fragmentation rate B(x) = B, which corresponds 
to the coefficient b(x,y) = 25 y=x / 2 : this was studied in [13], where exponential 
convergence of solutions to the equilibrium G was proved. Similar results were 
obtained for the mitosis case when B(x) is bounded above and below between 
two positive constants (i.e., for b{x, y) = 2B(x)S y=x / 2 )- An exponential speed of 
convergence has also been proved in [6] allowing for quite general fragmentation 
coefficients b, provided that the total fragmentation rate B{x) is a constant. Of 
course, in the above results existence of the equilibrium G and the eigenfunction 
4> were also proved as a necessary step to study the asymptotic behavior of 
equation (1). The problem of existence of these profiles was studied in its own 
right in [3, 8], where results are obtained for general fragmentation coefficients 
b, without the restriction that the total fragmentation rate should be bounded. 

The self-similar fragmentation equation The fragmentation equation is 

dtft = Cf t , (20a) 
Ms) = fin{x) (x > 0), (20b) 

which models a set of clusters undergoing fragmentation reactions at a rate 
b{x,y). In the study of the asymptotic behavior of ft, one usually restricts 
attention to fragmentation kernels which are homogeneous of some degree 7 — 1, 
with 7 > 0: 

b(rx,ry) = r 1 ~ 1 b(x,y) (r > 0, x > y > 0), (21) 

so B(x) = Bqx 1 for some -Bo > 0. Then one looks for self-similar solutions, this 
is, solutions of the form 

f t {x) = {t+lfh G {{t + l) 1 hx), (22) 

for some nonnegative function G. Observe that fo — G in this case. We omit 
the case 7 = 0, as self-similar solutions have a different expression in this case, 
and the asymptotic behavior of the fragmentation equation is also different, and 
must be treated separately. Such a remarkable function G is called a self-similar 
profile and is solution to the eigenvalue problem (5) with a(x) — x and A = 1. 
In this case one can easily show that in fact A = 1 is the largest eigenvalue of 
the operator g 1— > — d x (x g) + Cg. Existence of solutions G of equation (5)in this 
setting has been studied in [4] and also in [3] . The corresponding dual equation 
(6) is explicitly solvable in this case, with <f>(x) = Cx for some normalization 
constant C . 
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The above suggests to define g through the following change of variables: 

f t (x) = (t+l) 2 ^g(^-log(t+l),(t + l) 1 ^x) (t,x>0), (23) 

or, writing g t in terms of /, 

gt{x) := e- 2t /(e 7 * - 1, e^x) (t, x > 0). (24) 

Then, g t satisfies the self-similar fragmentation equation: 

d t g t + x d x g t + 2g t = 7 Cg t (25a) 
go(x)=f in (x) (x>0). (25b) 

We may redefine b(x,y) to include the factor 7 in front of Cg y , and omit 7 
in the equation. Then, this equation is of the form (1) with a(x) = x and 
A = 1. The long time convergence (9) or (10) means here that generic solutions 
asymptotically behave like the self-similar solution (t + l) 2 / 7 G((i + l) 1 / 7 ^). 

The problem of convergence to self-similarity for the fragmentation equation 
(20) was studied in [4], and then in [10]. Results on existence of self-similar 
profiles and convergence of solutions to them, without a rate, are available in 
[4, Theorems 3.1 and 3.2] and [10, Theorems 3.1 and 3.2], and are obtained 
through the use of entropy methods. To our knowledge, no results on the rate 
of this convergence were previously known. 



Assumptions on the fragmentation coefficient We turn to the precise 
description of our results, for which we will need the following hypotheses. 

Hypothesis 1.1. For all x > 0, b(x, ■) is a nonnegative measure on the interval 
[0,a;]. Also, for all ip g Co([0, +00)), the function x H> J, Q ,b(x,y)ip(y) dy is 
measurable. 

Hypothesis 1.2. There exists k > 1 such that 

[ b(x, y) dy = kB{x) (x > 0). (26) 

Hypothesis 1.3. There existO < B rn < Bm and"/ > —1 (for growth-fragmentation) 
or 7 > (for self-similar fragmentation) such that 

B m x^ < B(x) < Bmx 1 {x > 0). (27) 

Hypothesis 1.4. There exist /x, C > such that for every e > 0, 

b(x,y)dy <Ce^B(x) (x > 0). 

Hypothesis 1.5. For every S > there exists e > such that 

[ b{x,y)dy <SB(x) (x > 0). (28) 

J{l-e)x 
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Both Hypothesis 1.4 and 1.5 are "uniform integrability" hypotheses: they 
say that the measure b(x, y) dy is not concentrated at y — or y = x, in some 
quantitative uniform way for all x > 0. For some parts of our results we will 
also need the following hypothesis. 

Hypothesis 1.6. The measure b(x, •) is (identified with) a function on [0,x], 
and there exists Pm > ( actually, it must be Pm > 2 due to the definition of 
B) such that 

b(x,y)<P M — ( x >y>0). (29) 
x 

Hypothesis 1.7. There exists P m > such that 

b{x,y)>P m ^- (x>y>0). (30) 

Wc list the above hypothesis in a rough order of least to most restrictive. 
When taken all together they can be summarized in an easier way: Hypotheses 
1.1—1.7 are equivalent to assuming Hypothesis 1.1, 1.2 and that there exists 
< B rn < E>m satisfying 

2B m x 7 " 1 < b(x, y) < 2B M x'^ 1 (0 < y < x) (31) 

for some 7 > —1, or 7 > Oin the case of growth- fragmentation. We give 
Hypotheses 1.1-1.7 instead of this simpler statement for later reference: some 
of the results to follow use only a subset of these hypotheses, and it is preferable 
not to use the more restrictive conditions where they are not needed. 

As an example of coefficients which satisfy all Hypotheses 1.1—1.7 we may 
take b of self-similar form, i.e., 

b(x,y) -x^pf^j (x>y>0), (32) 

where p : (0,1) — > (P m ,PM) is a function bounded above and below, and 7 
within the specified range (for the main theorem below one must have 7 £ (0, 2] 
for self-similar fragmentation, and 7 G (0,2) for growth- fragmentation) . A bit 
more generally, one can take 

b(x,y) :=h{x)x^- 1 p{^j (x>y>0), (33) 

with h(x) : (0, +00) — > [B m7 B]\i] a function bounded above and below. 

Previous results We recall the following result, readily deduced from [3, 
Theorem 1 and Lemma 1] and [4, Theorem 3.1]: 

Theorem 1.8. Assume Hypotheses 1.1-1.4- There exists a unique triple (A, G, <fi) 
with A > 0, G G L 1 and <f> G W^^° which satisfy (5) and (6) (in the sense of 
distributional solutions for (5a), and of a.e. equality for (6a)). 



8 



In addition, 



G(x) > 0, <f){x) > (x > 0), (34) 
for all a > 0, x h-> x a a(x)G(x) e W 1,x (0, +oo). (35) 

in £/ie case of </ie growth-fragmentation equation (a(x) = 1) it holds that 
0(0) > 0. 

In the case of the self- similar fragmentation equation (a{x) = x), assume in 
addition that, for some k < and C > 1, 

y fc 6(a;, y) dy < C x k B(x) (x > 0). (36) 

Jo 

TTien, 

/>0O 

x k G(x) dx < +oo. (37) 



The above result can actually be proved under weaker hypotheses on the 
fragmentation coefficient b (see [3]) but we will only use the given version. We 
note that (37) is derived in [4] under more restrictive conditions on b (in partic- 
ular, it is asked there that b is of the self-similar form (32)); however, the same 
proof can be followed just by using (36), and we omit the details here. 



Main results The following theorem gathers our main results in this work. 

Theorem 1.9. Consider equation (1) in the self-similar fragmentation case 
(a(x) = x) or the growth-fragmentation case (a(x) = 1), and assume Hypotheses 
1.1-1.7. Also, 

• For self- similar fragmentation, assume 7 £ (0,2]. 

• For growth-fragmentation, assume 7 £ (0, 2) or alternatively that for some 
B b > , 

b(x,y) = ^ (x>y>0). (38) 
x 

Denote by G > the asymptotic profile (self-similar profile in the first case, first 
eigenfunction in the second case) as well as by <f> the first dual eigenfunction 
(4>{x) = x in the first case). These equations have a spectral gap in the space 
L 2 (4> G~ l dx). 

1. More precisely, there exists f3 > such that 

VgeX, H 2 (g\G)<^D\g\G), 
(the right hand term being finite or not), where we have defined 

H := L 2 (G~ X (j>), X:={g£H; J ^ = ^ G0=i}. 
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Consequently for any gi n £ X the solution g £ C([0, oo); Lh) to equation 
(1) satisfies 

Vt>0 \\g t -G\\ H <e- 2 ^\\g in -G\\ H . 

2. For this second part, in the case of growth-fragmentation, we need to as- 
sume additionally that 7 £ (0,1) (no additional assumption is needed for 
the self-similar fragmentation case). There exists k = £(7, Pm) > 3 and 
for any a £ (0,/3) and any k > k there exists C a ,k > 1 such that for any 
gin £ % := L 2 {9), 9{x) = <j>(x) + x k , J 4>g m = 1, there holds: 

Vt > \\g t -G\\ H < C a>k e~ at \\g in - G\\ H . 

The main improvement of this result is that it allows us to treat total frag- 
mentation rates B which are not bounded, for which no results were available. 
Also, the result is obtained through an entropy-entropy dissipation inequality, 
which is not strongly tied to the particular form of the equation, and may be use- 
ful in related equations. We also observe that this is to our knowledge the first 
result on rate of convergence to self-similarity for the fragmentation equation. 

In our result we include just one case where B is a constant, with b given 
by the particular form (38) in the case of growth-fragmentation. The difficulty 
in allowing for more general b (with B still constant) is not in the entropy- 
entropy dissipation inequality, but in the estimates for the solution <fi of the 
dual eigenproblem (6). The estimates in section 4.2 are not valid in this case, 
and in fact one expects that <f> should be bounded above and below between 
two positive constants. However, we are unable to prove this at the moment, 
and anyway the estimates of <f> are not the main aim of the present paper. 
Consequently, we state the result only for the b in (38), for which <f> is explicitly 
given by a constant. (We notice that for self-similar fragmentation, the estimates 
on G fail for 7 = 0, so this case is not included). 

On the other hand, the positivity condition (30) is strong, and makes this 
result not applicable to some cases of interest, such as the mitosis case mentioned 
at the end of section 1. Notice that, as remarked at the end of section 2 in [6], 
entropy-entropy dissipation inequalities are actually false in this case, which 
shows that a different method is required for this and similar situations where 
the fragmentation coefficient is "sparse" enough for the inequalities to fail. 

Strategy The main difficulty for establishing Theorem 1.9 lies in proving point 
1, while point 2 is a consequence of point 1 and of a method for enlarging the 
functional space of decay estimates on semigroups from a "small" Hilbert space 
to a larger one recently obtained in [5] . 

The strategy we follow to prove inequality (17) is inspired by a paper of 
Diaconis and Stroock [2] which describes a technique to find bounds of the 
spectral gap of a finite Markov process. Equation (1) can be interpreted as the 
evolution of the probability distribution of an underlying continuous Markov 
process, and as such can be thought of as a limit of finite Markov processes. 
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The expression of the entropy and entropy dissipation has many similarities to 
the finite case, and hence the ideas used in [2] can be extended to our setting. 
The basic idea is that one may improve an inequality in the finite case by 
appropriately choosing, for each two points i and j in the Markov process, 
a chain of reactions from i to j which has a large probability of happening, 
measured in a particular way. In our case this corresponds to choosing a chain 
of reactions that can give particles of size x from particles of size y, and which 
has a better probability than just particles of size y fragmenting at rate b(y, x) 
to give particles of size x. In fact, we need the additional observation that 
one can choose many different paths from y to x, and then average among 
them to improve the probability. These heuristic ideas are made precise in the 
proof of point 1 in Theorem 1.9, given in section 2. This suggests that the 
same techniques may be applicable to other linear evolution equations which 
can interpreted as the probability distribution of a Markov process. 

In the context of the Boltzmann equation, an idea with some common points 
with this one was introduced in [1] to prove a spectral gap for the the linearized 
Boltzmann operator. However, the inequality needed in that case does not 
have the same structure since the linearized Boltzmann operator does not come 
from a Markov process (it does not conserve positivity, for instance), and the 
"reactions" to be taken into account there are not jumps from one point to 
another, but collisions between two particles at given velocities, to give two 
particles at different velocities. Nevertheless, the geometric idea does bear some 
resemblance, and the connection to techniques developed for the study of finite 
Markov processes seems interesting. 

Concerning the proof of point 2, we mainly have to show that the operator 
involved in the fragmentation equation decomposes as A + B where A is a 
bounded operator and B is a coercive operator in the large Hilbert space T-L. 
Then, for a such a kind of operator, the result obtained in [5] ensures that the 
operator in the large space T-L inherits the spectral properties and the decay 
estimates of the associated semigroup which are true in the small space H. 
The additional restriction 7 £ (0, 1) in the extension of the spectral gap comes 
from the fact that we were unable to prove the coercivity of the operator B for 
7 6 [1,2). 

In the next section we prove the functional inequality (17). In sections 3 and 
4 we prove upper and lower bounds on the profiles G and 4>, solution of the eigen- 
problcm (5)-(6); section 3 is dedicated to the self-similar fragmentation equation 
(the case a(x) = x), while section 4 deals with the growth- fragmentation equa- 
tion (the case a(x) = 1). Finally, section 5 uses these results and the techniques 
in [5] to complete the proof of Theorem 1.9. 

2 Entropy dissipation inequalities 

In order to study the speed of convergence to equilibrium, in terms of the en- 
tropy functional, we are interested in proving the entropy-entropy dissipation 
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inequality (17): 



H 2 [g\G}<±D»[g\G}. 



We begin with the following basic identity: 

Lemma 2.1. Take nonnegative measurable functions G,cj> : (0, oo) — > R + such 
that 

pOO 

/ <p(x)G(x)dx = 1. (39) 
Jo 

Defining 

poo poo 

D 2 [g\G] := / / 0(s) G(s) G{y) (u(x) - u{y)) 2 dydx (40) 

JO Jx 

with u(x) = q(x) ' there holds 

H 2 [g\G] = D 2 [g\G] (41) 
for any nonnegative measurable function g : (0, oo) — > R such that J <fi g = 1. 

Proof. With the above notations, by a simple expansion of the squares, we find 
that, 

/•OO POO 

H 2 [g\G] = / (u(x) - l) 2 G(x) 4>{x) dx = / u{x) 2 G{x) 4>{x) dx - 1, 
Jo JO 

while expanding £>2[g|G] we find 



OO POO 



D2[g\G] = - I I G(x)cf>(x)G(y)4>(y)(u(x)-u(y)y dydx 



JO 

2 Jo 



oo 

2/ 



which gives the same result. □ 

Taking into account the above result, in order to prove (17) it is enough to 
prove the following inequality for some constant C > 0: 

D 2 [g\G] <CD b [g\G}. (42) 

Of course, one gets (42) if one assumes that the function appearing below the 
integral signs in the functional D 2 [g\G] is pointwise bounded by the corre- 
sponding function in the functional D b [g\G], or more precisely if one assumes 
G(x) 4>(y) < C b(y, x) whenever < x < y. However, the range of admissible 
rates 6 for which such an inequality holds seems to be very narrow. For example, 
for the self-similar rate (32), one must impose 7 = 2. 

The cornerstone of the proof of the functional inequality (42) lies in splitting 
D 2 [g\G] in two terms: one for which the pointwise comparison holds and one 
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where it does not. In the latter part, called Z?2,2[5|G] in the next lemma, the 
idea is to break u(x) — u(y) into "intermediate reactions" in order to obtain a 
new expression of D 2 [g\G] for which the pointwisc estimate applies. The main 
part of the proof of inequality (42) is then the bound for this second part of 
D 2 [g\G]. And it is just the content of the following result: 

Lemma 2.2. Take measurable functions (f>, G : (0, +00) — > R+ with 

G(x)4>(x)dx = l (43) 



and such that, for some constants K, M > and R > 1 and some function 
C : (R,oo) -> [l,oo), 

< G{x) <K (x> 0), (44) 
G(y)(f>(y) dy < KG(x) (x > M), (45) 

Rx 

-^-4>{y) < K 4>{z) (m&x{2RM,Rz} <y <2Rz). (46) 

Civ) 

Defining 

/>oo />oo 

D2.2 [g\G] := / / 4[x) G(x) <P(y) G{y) (u(x) - u{y)) 2 dy dx (47) 

JO J-aiax{2Rx.2RM} 

there is some constant C > such that 

POD pOO 

D 2i2 [g\G]<C / ay)y- 1 H^G(y)(u(x)-u(y)) 2 dydx 1 (48) 

JO Jiaax{x,M} 

for all measurable functions g : (0, 00) — )• R + (in the sense that if the right hand 
side is finite, then the left hand side is also, and the inequality holds). 

The point of this lemma is that the right hand side of (48) will be easily 
bounded by D b [g\G}. 

Remark 2.3. Inequality (48) is strongest when ((x) = 1. Moreover, we expect 
(46) to be true when £ is constantly 1 (for <j> the solution of the dual equation 
(6), under some additional conditions). The reason that we allow for a choice 
of £ is that we are not able to prove that (46) holds with £ = 1, but only for a 
function Q = Q(x) which grows like a power x £ , with e as small as desired. 

Proof. We will denote by C any constant which depends on G, <fi, K, M, or R, 
but not on g. 

First step. The idea is to break u(x) — u(y) into "intermediate reactions" : for 
y > x and any z € [x,y], one obviously has 

u(x) — u{y) = u(x) — u(z) + u{z) — u(y). 
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We then average among a range of possible splittings. More precisely, for y > 
2Rx, 



2R 



ry/R 

u(x) — u(y) = — / (u(x) — u(z) + u(z) — u(y)) dz, 

V Jy/(2R) 



and then by Cauchy-Schwarz's inequality, 
2R f y l R 

(u(x) — u(y)) 2 < — / (u(x) — u(z) + u(z) — u(y)) 2 dz 

V Jy/(2R) 

4R f v/R / / \ ,„ 2 , 4i? f y/R . . . . .. 2 , 

< — / (u(x) - u(z)) 2 dz -] / (u(z) - u(y)f dz 

=■ Ti{x,y) + T 2 {y). 

Using this, we have 

/>oo />oo 

D2.2 [g\G] < / / <f>(x)<t>(y)G(x)G(y) T x (x, y) dy dx 

JO Jmax{2Rx,2RM} 
roo poo 

/ <j ) (x)cj ) (y)G(x)G(y)T 2 (y)dydx. 

IQ Jmax{2Rx,2RM} 

We estimate these terms in the next two steps. 
Second step. For the term with Ti, 

00 POO 

/ 4>{x)4>(y)G{x)G{y)T r (x,y)dydx 

Jm&x{2Rx,2RM} 

poo poo 1 ry/R 

= C / cj){x)(j)(y)G{x)G{y)- (u(x) - u{z)f dz dy dx 

J0 Jma,x{2Rx,2RM} V J y / (2R) 

= C [ [ 4>{x)G{x){u{x) -u(z)) 2 [ ^-G{y)dydzdx 

JO Jmax{x,M} Jmax{2Rx,Rz,2RM} V 

pOO pOO -1 />OG 

<C / 4>{x)G{x){u{x)-u{z)) 2 — / cj ) (y)G(y)dydzdx 

JO Jmax{x,M} Kz J Rz 

/• 00 /• 00 1 

<C / <f>(x)(u(x) - u(z)) 2 -G{z)dzdx 

JO Jma,x{x,M} z 

where we have used (45) and the fact that G is uniformly bounded from (44). 
With this series of inequalities we obtain that the last one integral is bounded 
by the right hand side in (48), as £(y) > 1 for all y. 

Third step. For the term with T 2 , we use (43) and (46). Note that (46) was 
not used before this point in the proof: £(y) < 1 was enough up to now, but 
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this is not the case in the following calculation: 

<j>(x) ( / ) (y)G(x)G(y)T 2 (y)dydx 



OO />OC 



Jmax{2Rx,2RM} 

oo poo -| py / R- 

' ' ' - ' ^2 



= C / 4>{x)4>(y)G{x)G{y)~ (u(z)-u(y)Ydzdydx 

Jo Jmax{2Rx,2RM} V Jy/(2R) 

roo r y/R i r y/(2R) 

= C / -<b(y)G(y)(u(z) - u(y)) 2 / #c)G(z) dj, 

./2.RM Jy/(2R) V JO 
poo py/R i 

<C / / -4>(y)G(y)(u(z) - u{y)) 2 dzdy 

J2RM Jy/(2R) V 

= C r ^ £M G (y) («(*) - u(y)) 2 dy dz 

JM Jmax{2RM,Rz} V ^(U) 

<c r f z M^ z)G{y){u{z) _ u{v) f dydZj 

Jo Jmax{2RM,Rz} V 



where the last integral is bounded by the right hand side in (48) and this finishes 
the proof. □ 

Once we have controlled the "bad" term in D 2 [g\G] we can reach the ob- 
jective of this section: to obtain an entropy-entropy dissipation inequality. It is 
shown in the following theorem: 

Theorem 2.4. Assume the conditions in Lemma 2.2, and also that the frag- 
mentation coefficient b satisfies, for the constants K,M,R in Lemma 2.2, 

G{x) 4>{y) < K b(y, x) (0 < x < y < max{2i?x, 2RM}), (49) 
C(y) ZT 1 < K b(y, x) (y > M, y>x> 0). (50) 

Then there is some constant C > such that 

H 2 [g\G]<CD b [g\G] (51) 

for all measurable functions g : (0, oo) — > R+ such that f^°g<ft = 1 (in the 
sense that if the right hand side is finite, then the left hand side is also, and the 
inequality holds). 

Proof. We split D 2 [g\G] as 

D 2 [g\G] = D 2A [g\G]+D 2a [g\G], 

with 

poo pnmx{2Rx,2RM} 

D2.1 [g\G] := / / 4>(x) G{x) 4>{y) G(y) (u(x) - u{y)f dydx 

JO Jx 
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and £>2,2 [g\G] denned by (47). On the one hand thanks to (49) we have 



D 2 ,i i9\G]<j o J K b(y, x)<f>(x) G(y) (u(x) - u{y)) 2 dydx 

< KD b [g\G}. 

On the other hand, thanks to inequality (48) in Lemma 2.2 and (50) we have 

p OO pOO J-f \ 

D 2 ,2[g\G}) < C / ^-c/>(x)G(y)(u(x)-u(y)) 2 dydx 

JO Jmax{x,M} V 

p OO pOO 

< CK / b{y,x)<t>(x)G(y)(u(x)-u(y)) 2 dydx 

JO Jmax{x,M} 

< CKD b [g\G}. (53) 

We conclude by gathering (52) and (53). □ 

This result will be the key to prove point 1 in Theorem 1.9, which will be 
done in section 5.1. In order to prove that (44) and (49) hold in our context, 
we need to establish some upper bounds on G and (j) for our model cases. For 
(45) and (46) we need to control the asymptotic behavior (bounds from above 
and below) of the functions G(x) and 4>{x) as x — > oo. These upper and lower 
estimates on G and <p will be established in the next sections. Finally, condition 
(50) simply imposes some restrictions on the fragmentation rate (typically some 
restrictions on the value of 7 for a self-similar fragmentation rate of the form 
(32)). 



3 Bounds for the self-similar fragmentation equa- 
tion 

In order to apply Theorem 2.4 to the self-similar fragmentation equation (25) 
we need more precise bounds than those proved in [4, 10, 3]; in particular, we 
need L°° bounds on the self-similar profile G for condition (44) to hold. Wc 
actually prove the following accurate exponential growth estimate on the profile 
G. 

Theorem 3.1. Assume Hypotheses 1.1-1.5 on the fragmentation coefficient b. 
CallA{x) :=J*^ds. 

1. For any S > and any a G (0, Bm/Bu), a' G (1, +00) there exist constants 
C' = C'(a', 6), C = C(a, S) > such that 

C' e r a ' A( - x) < G(x) < C e"° A W forx>5. (54) 

2. Assume additionally Hypothesis 1.6. Then one may take C independent 
of 5 in (54 ), i.e., 

G(x) < Ce' aA{x) for x > 0. (55) 



16 



Remark 3.2. We notice that, due to Hypothesis 1.3, (£> m /7)a; 7 < A(x) < 
(Bm/^x 1 , so the bound (54) directly implies that for any a\ > Bm/j and 
a 2 < Bf n /(^B M ) one has 

C e~ ai xl < G{x) < C e~ a2 x ~* for x > 6, (56) 

and also the corresponding one instead of (55). Note that when Bm = B m the 
condition on oi, a 2 becomes a 2 < BmIi < Oi- 

The goal of this section is to give the proof of the above theorem, which we 
develop in several steps. We remark that in the particular case of p constant 
in (32) (so, p = 2 due to the normalization J zp(z) dz = 1) we can find an 
explicit expression for the equilibrium G (i.e., a solution of (5)). Indeed, G{x) = 

cxp ^— S^il ds^J satisfies (5). As a consequence, in the case b(x,y) = 2x 1 ~ 1 

(so B{x) = x 7 ), the profile G is e~~ for 7 > 0. In the general case where b is 
not of the form (32) there is no explicit expression available for the self-similar 
profile G. 

Lemma 3.3. Assume that b satisfies Hypotheses 1.1-1.5, and consider G the 
unique self-similar profile given by Theorem 1.8. For any < a < B m j^ there 
exists a constant C a such that: 



f 

Jo 



e ax '' G(x) dx < C a . (57) 



Proof. We denote by Mfc the fc-th moment of G. Multiply eq. (5) by x k with 
k > 1 to obtain: 



00 

k ; 



(k - l)M k > (1 - Pk ) / x k B(x)G(x) dx. (58) 
Jo 

where we have taken into account that, using Corollary 6.4, 

1 /*oo /*oo py 

x k b(y,x)G(y)dydx = G(y) x k b(y,x)dxdy 

Jx JO Jo 

poo 

< Pk y k B(y) G(y) dy. 
Jo 

By Hypothesis 1.3, and as pk < 1 for k > 1, (58) implies that 

k 

(k - 1) M k > (1 - p k )B m M k+7 , or M k+1 < — — M k . 

B m {l -pk) 

Applying this for integer £ > 1 and k := 1 + £j, 

1 + 7^ 



M 1+( ^ +1)7 < CiM 1+ill where C e 



B m {l —Pl+yt) 
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Solving the recurrence relation, 

e-i 

M l+h < M 1+1 JJ d (£>1). (59) 

i=l 

With this, 

/ e^Gix) dx = V / — x 7J G(a;) = V — M 7 , 
l= o ^ z - i=o *" 

< M + M 7 + M 1+7 £ ^ [] Ci. 



The last expression in the sum is a power series in a, with radius of convergence 
equal to B m /j. This can be checked, for example, by noticing that 

as t: — > +oo, 



£ + 1 B„ 

which corresponds to the quotient of two consecutive terms in the power series. 
This proves the lemma. □ 

With this result we can now prove Theorem 3.1. 
Proof. Taking 8 > 0, we give the proof in several steps. 

First step (upper bound for x > 5): From (5), we calculate as follows: 
d x (x 2 e aA G) = x e aA (2 G + xd x G) + a B x e aA G = xe aA CG + aBx e aA G 

p OO 

= x e aA I b{y, x) G{y) dy + (a - 1) B x e aA G. (60) 

J X 

Let us show that for a < 1 this expression is intcgrablc. For the first term, 

i poo poo py 

xe aA{x) / b{y 7 x)G{y)dydx = / G{y) / e aA( - x h{y , x) dx dy 

Jx JO JO 

< y e aA ^G(y)B(y) dy < B M / y^ +1 G(y) dy < +00, 

Jo Jo 

where the last inequality is due to Lemma 3.3, since y 1+1 e 1 ^ < C e by ~' , 
where a < b < and C is a constant depending on a, b and 7 (recall that a < 
B m /BM)- For the same reason, the last term in (60) is integrable. Therefore, 
since d x (x 2 e aA ^G(x)) is bounded in L 1 , we deduce that 

x 2 e axl G{x) E W(0,oo) C L°° 

and in consequence 

G(x) < Cler aA ^ for x > 5. 
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Second step (lower bound for x > 5): Writing equation (60) for a = 1 
gives that x \-> x 2 e A ^G(x) is a nondecreasing function, so 

x 2 e A{x) G{x) > 5 2 e A ^G{6) (x > 6), 

which implies the lower bound in (54) for any a < 1. Notice that G(6) > by 
Theorem 1.8. 

Third step (upper bound for x < 5): As x n> x 1 e x ~' 1 1 G{x) has bounded 
variation, it must have a limit as x — > 0. As x i— » x G{x) is intcgrablc, it follows 
that this limit must be 0. Hence, writing (60) for a = 1, integrating between 
and z, and using Hypothesis 1.6, 

/>Z poo 

z 2 e z ~' h G{z) = xe aA b(y 7 x) G(y) dy dx 

JO Jx 

f-Z f-OO f-Z 

< CB M / xe aA / y-*- 1 G(y) dydx < C xe x ~ ' h dx < C'z 2 e z ~'^, 

JO Jx Jo 

which implies G(z) < C for all z > 0. This is enough to have (55) for x < 
S. We have used above that the moment of G of order 7 — 1 is bounded, 
as given by Theorem 1.8, taking into account that (29) holds and 7 > 0, so 
f*y-y- 1 b(x,y)dy<(C/ 1 )x'r-\ □ 



4 Bounds for the growth-fragmentation equa- 
tion 

In this section we present some estimates by above and below for the functions 
G and 4>, solutions to the eigenvalue problem (5) and (6). The aim, as in 
the previous section, is to obtain bounds which are accurate enough to apply 
Theorem 2.4, and then prove Theorem 1.9. 

The main additional difficulty as compared to the self-similar fragmentation 
equation from previous section is that the dual eigenfunction 4> is in general not 
explicit, which makes it necessary to have additional estimates for it. 

In the rest of this section we will give the proof of the following result: 

Theorem 4.1. Assume Hypotheses 1.1-1.5. Call 

K(x) :=\x+ [ B(x) dx. 
Jo 

1. For any a S (0, B m / Bm) there exists K a > such that 

Vx > G(x) < K a e - aA W. (61) 

If we also assume Hypothesis 1.6 and"/ > 0, then this can be strengthened 
to 

Vx>0 G(x) <K a mm{l,x}e- aA{x) . (62) 
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2. For any 5 > there exists K$ > such that 

Vx>5 G(x)>K s e- A( - x) . 



(63) 



3. Assume additionally that 
B(x] 



x^ 



and B(x) — > +oo as x — > +oo. (64) 



(Here fi is the one in Hypothesis 1-4-) There exist Co,Ci £ (0, oo) and 
for any k 6 (0, 1) there exists C k £ (0, oo) such that 

Vx>0 C + C k x k < cf>{x) < Ci (1 + x). (65) 

Remark 4.2. With (61), (63) and (27), in the case 7 > it is easy to see that for 
any 5 > 0, ai > B M /{l + 1) and a 2 < B 2 m l(B M (l + 1)) there exist Ci, C 2 > 
such that 

Ci e" ai " 7+1 < G(z) < C 2 e" a2 (x > S). (66) 

In the case 7 = 0, from (61) and (63) one sees that for any S > and a > 1 
there exist C\ , C 2 > such that 

Cie -(\+B M )x < G(a;) < C r ae -.(A+fl m )x (x > (g 7) 

4.1 Bounds for G 

Lemma 4.3. Assume hypotheses 1.1-1.5. Call 

A m (x) := Ax+-^^ +1 . 

7+1 

For eac/i a < 1, 

/■OO 

e oAm ^)G(x)dx < +00. (68) 

Proof. We will first do the proof for 7 > 0, and leave the case — 1 < 7 < for 
later. Multiply equation (5) by x k with k > 1, and integrate to obtain 

/•oo />y poo 

kM k _ 1 -\M k + G(y) x k b{y,x)dxdy- B(x) x k G dx = 0, 
Jo Jo Jo 

which gives, using (92) and then (27) (noting pf. < 1 for k > 1), 

(1 - p k )B m M 1+k < feM fc _! - AMfc < fcM fc -i. (69) 
Applying this for k = £(j + 1) — 7, with £ > 2 an integer, 

ih + 1) - 7 

M eh+1) < CfAf (w)(7+ i), with C e := — * f— . 

(1 -p^( 7+1 )_ 7 J±f m 
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Solving this recurrence relation gives, for I > 2, 

t 

M £(7+1) < M 7+1 Y[C t (£>2). 

i=2 

Now, following an analogous calculation to the one in Lemma 3.3, 

/■oo 00 i £ 

/ e a G{x) dx < M + M 7+1 ^ — ]J C e . (70) 

J ° 1=2 l - t=2 

Again as in Lemma 3.3, one can check that the above power series in a has 
radius of convergence B m /(-f + 1) (using that pk — > when k — > +00, from 
Corollary 6.4 in the Appendix). This proves the lemma for 7 > (note that the 
dominant term in A in this case is a; 7+1 , as then x < C e + ex 1+1 for any e > 
and some C c > 0). 

When 7 = 0, from (69) we obtain 

((1 - Pk )B m + X)M k < kM k - X (k > 1). 

We can then follow the same reasoning as above, with the only difference that 
now 

C '^ (l-P,) B ... + A <^ 2 '- 

Now the power series in (70) has radius of convergence B m + A, which proves 
the lemma also in this case. 

In the case 7 G (—1,0), from (69) we obtain the inequality 

Vfc>l M k <jM k _ u 

from which we deduce thanks to an iterative argument as before that 

V/ceN* M k <-^(\Mi). 

Hence, for a < 1, 

/ e aXx G(x) dx = M Q + ^ -rp A/ fe < M o + AAfi ^ a k < +00. 

J® k=l ' k=l 

The dominant term in A in this case is Ax, so this finishes the proof. □ 



Proof of points 1-2 in Theorem J^.l. With the previous lemma we are ready to 
prove our bounds on G. 
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First step (upper bound): Take < a < 1, and calculate the derivative of 

G(x)e aA ^: 

(Ge aA Y = (a - 1)(B + X)Ge aA + e aA C+{G). (71) 

For a < 1 one can see that the right hand side is integrable on (0, +oo): for the 
last term, 



e aA C+{G)= / C* + (e aA )G= / G(x) / e aAfe; ] b{x , y) dy dx 
Jo Jo Jo 

/>oo f-X f-OO 

< G{x)e aA(x) b(x,y)dydx = K B(x)G(x)e aA ^ dx < +oo, 



where we have used (26). The last expression is finite for a < B m /BM due to 
Lemma 4.3 and the fact that 



A(x) = Xx+ / B(y)dy < Xx + B M 

Jo Jo 



V 1 dy 



, . B M +1 B M f, B m +1 \ B M . / v 
= Xx-\ — r^ 7 <— AaH — -x 7+ = — — A m (x), 

7+1 B m \ 7+1 / B m 

using the upper bound in (27). The other term in (71) is also integrable for 
similar reasons, and we deduce that e aA G e W(0,oo) C L°°, which proves 
(61). 

In order to get (62) we need to prove that additionally, G(x) < Cx for x 
small (say, x < 1). For this it is enough to notice that G(0) = due to the 
boundary condition (5b), and also that the right hand side of (71) is bounded 
for x £ (0, 1), as B and e aA G are, and 

C+(G) = / b(y, x)G(y) dy < C y^ l G{y) dy < +oo, 

Jx JO 

due to Hypotheses 1.6 and 1.3. The last integral is finite because 7 > and we 
already know G is bounded. This finishes the proof. 

Second step (lower bound): Writing equation (71) for a = 1 we obtain 

(Ge A )' = e A £+(G) > 0, 

and hence 

G{x)e A{x ^ > G(S)e A ^ (x > 6). 
This proves the result, as G(6) > by Theorem 1.8. □ 

4.2 Bounds for <fi 

First, in the case B(x) = B constant, the first eigenvalue A of the operator 
— d x + C is explicitly given by A = B(k — 1) (under Hypotheses 1.1-1.2), and 
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4>{x) = C constant is a solution of (6), for some appropriate C > determined 
by the normalization (6b). 

In general the solution <f> of the eigcnproblcm (6) is not explicit, and for its 
study we will use the following truncated problem: given L > 0, consider 

-d x cf> L + (B(x)+\ L )(t> L (x) = £* + (4> L )(x) (0<x<L, (72a) 
(f>>0, (f> L (L) = 0, [ G(x)<j> L {x)dx = 1. (72b) 

This approximated problem is slightly different from the one considered in [3], 
in that we arc considering the first eigenvector G in the normalization (72b), 
and not an approximation Gl obtained by solving a similar truncated version 
of equation (5). However, this modification is not essential, and the results in 
[3] show the following (see also the truncated problems in [8], [12]): 

Lemma 4.4. Assume Hypotheses 1.1-1.4- There exists Lq > such that for 
each L > Lq the problem (72) has a unique solution (Xl,<Pl), with Xl > and 
4>L £ Wl^ . In addition, 

X L L ^°° A, (73) 
for every A > 0, 4>l L zHf° ^ uniformly on [0, A), (74) 

where {X,(j)) is the unique solution of (6). 

In the rest of this section we always consider L > Lq, so that Lemma 4.4 
ensures the existence of a solution. 



4.2.1 Upper bounds 

In order to obtain bounds for (j> we use a comparison argument, valid for each 
truncated problem on [0, L], and then pass to the limit, as the bounds we obtain 
are independent of L. Let us do this. The function <pL is a solution of the 
equation 

Scj> L {x)=Q (xe(0,L)), 
where S is the operator given by 

S<j>(x) := -cj>'{x) + (X L + B(x)) <j>(x) - f b(x, y) <j>{y) dy, (75) 



dchned for all £ W 1,oo (0, L), and for x £ (0, L). The operator S satisfies the 
following maximum principle: 

Definition 4.5. We say that w £ W 1,00 (0, L) is a supersolution of S on the 
interval / C (0, L) when 

Sw(x) > {x el). 
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Lemma 4.6 (Maximum principle for S). Assume Hypotheses 1.1. Take A > 
1/Xl- If w is a supersolution of S on (A,L), w > on [0,A] and w(L) > 
then w > on [A, L] . 

Proof. We will prove the lemma when w G C 1 ([0, L]), and then one can prove it 
for w £ W 1,oc (0, L) by a usual approximation argument. Assume the contrary: 
there exists xo £ (A,L) such that w(xq) < and w(x)/x attains a minimum, 
i.e., 

w{x ) < 0, (76) 

!^2<^ (xe(0,L)). (77) 

a; 

Then, because of (77), we have w'(xo) = w(xq)/xq and hence 

w(xo) r° 

S(x ) = h (Al +B(x )) w(x q ) - / b(x ,y)w{y)dy 

Jo 

W(X ) w(x ) f X ° 

< h (Al + B{xq)) w(xo) / b(xo,y)ydy 

Xo Xo Jo 



w(x )(\ L - — ) <0, 

V Xn) 



Xo> 

which contradicts that w is a supersolution on [A, L). □ 

One can easily check that v(x) = x is a supersolution of S on (1/Al, £). A 
useful variant of that fact is the following: 

Lemma 4.7. Assume Hypothesis 1.1-1.4, an d a ^ s0 that 

B(x) 

as x -» +oo, (78) 

where \x is that in Hypothesis 1.4- Take a smooth function 77 : [0,+oo) — > [0, 1], 
with compact support contained on [0, R] . Then, there exists A > and > A 
such that the function 

v(x) = x + r](x) 
is a supersolution of S on [A, L) for any L > L*. 
Proof. We calculate, using (4), 

Sv(x) = -1 + \ L x - rf[x) + (Al + B(x))r](x) - / b(x, y)rj(y) dy, 

Jo 

which for x > R becomes 

Sv(x) = -1 + X L x - J b(x, y)r,(y) dy > -1 + X L x - CR^^ 

as 77 is bounded by 1, and using Hypothesis 1.4. Due to (78), this is positive for 
all x greater than a certain number A which depends only on Al- As Al — > X 
when L —> +00 (see Lemma 4.4), one can choose A to be independent of L. □ 
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Proposition 4.8. Assume the hypotheses of Lemma J^.l. The solution <j> of (6) 
satisfies 

4>{x)<C{l + x) (x>0) (79) 

for some C > 0. 

Proof. Due to the uniform convergence (74) we have the bound 

(j)(x)<K(A) (ie[0,A]) (80) 

for some constant K = K(A) which does not depend on L. This proves the 
bound on (0, A] for a fixed A. 

To prove the bound on all of (0, +oo) consider the function v(x) = x + w{x) 
from Lemma 4.7 with n = l[o,i], which is a supersolution on (A,L) for some 
A > 0. Then, for any C > 0, the function 

w(x) := Cv(x) — 4>(x) 

is a supersolution on (A, L). Since <fr is bounded above on (0,A) by (80), uni- 
formly in L, we may choose C > ||<^||l=°(o,a) independently of L, such that 

<jj(x)<Cv(x) (xe[0,A]), 

or equivalently w > on [0,A]. As 4>(L) = 0, so that w(L) > 0, Lemma 4.6 
shows that to > on [0, L], which is the bound we wanted. □ 

4.2.2 Lower bounds 

Let us look now for subsolutions. 

Lemma 4.9. Assume Hypothesis 1,1—1,5, and also that B(x) — > +00 as x — > 
+00. Let if> : (— oo,0) —> [0,1] be a decreasing C 1 function which is 1 on 
(— 00, — e), on (— e/2,0), and satisfies \ip'(x)\ < 4/e for x G (— 00, 0). Take 
< k < 1. There is a number A which is independent of L for which v(x) := 
x k (p(x — L) is a subsolution of S on (A, L). 

Proof. First, from Corollary 6.4 we have 




y k b(x,y)dy>p' k x k B(x) (x > 0) 



for some p' k > 1. Hence, by Lemma 6.3, we may choose e» > such that 



f (1 e ' )X y k b(x,y)dy>C k x k B(x) (x>l), (81) 
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with Cfc > 1. Now, define ipl(%) = f{ x — L). We have, for x > max{i?i,, e/e»}, 
and using (81), 



Sv(x):=-kx k i <p L (x)-x k (p' L (x) + (\ L +B(x))x k <p L (x)- b(x,y) y k <p L (x) dy 



o 



<-x k + (X L + B{x))x k - I b(x, y) y k dy 



x k (- + \ L -B{x)(C k -l)), 



which is negative for x greater than some number A which depends only on Ai 
and k. In order to be able to apply (81) we have also used that x — e > (1 — e*)x 
for x > e/e*. □ 

Lemma 4.10. Assume Hypothesis 1.1—1.5, and also that B(x) — > +oo as x — > 
+oo. For any < k < 1 there is some constant Ck > such that 

4>(x) > C k x k (x > 0). 

Proof. Take if as in Lemma 4.9, and let A be the one given there. The function 

w(x) = 4>l{x) - Cx k ip(x - L) 

is a supersolution on (A, L) for any choice of C > 0. Now, the uniform conver- 
gence of {4>l} from Lemma 4.4 together with the positivity of </> from Theorem 
1.8 imply that there exists Ca such that for L large enough 

M*)>Ca (xG[0,A}), 

which in turn implies 

H{x) > (A- k C A )x k (^[0,4]). 

With C := A~ k Ca we have w > on [0, A], w(L) = 0, and we conclude 
using the maximum principle from Lemma 4.6, and again the locally uniform 
convergence of {4>l} from Lemma 4.4. □ 

5 Proof of the main theorem 

Finally, we are ready to complete the proof of Theorem 1.9. We give the proof 
of point 1 in the next subsection, and that of point 2 in the following one. 

5.1 Exponential convergence to the asymptotic profile 

Now that the inequality in Theorem 2.4 and the bounds on the profiles G and 
(p have been shown, we can use them to prove the first point in Theorem 1.9 
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Proof for self-similar fragmentation. Let us show that equations (44) , (45) , (46) , 
(49) and (50) hold for G, <j>. Then, as a direct application of Theorem 2.4, point 
1 of Theorem 1.9 follows. 

• The bound (44) is a immediate consequence of (56). With ((y) = 1 and 
whatever M and R be, (46) is satisfied due to the fact that cj>(y) = y for 
the self-similar fragmentation model. 

• Using (56) for any a\ > Bm/i and ai < B^/(7 Bm) (with < 7 < 2) we 
have for any x > M, R > 1 and for some constants denoted by C 

/>oo />oo />oo 

/ yG(y)dy<C ye^^dy^C y^ 1 e~ a V dy 

JRx JRx JRx 

= Ce- a ' Rlx ~' < Ce- aix '' < CG{x), 

where we consider a 1 such that: < a' < a-2 (< ai) (which is possible 
since R > 1), and that proves (45). 

• We split the proof of (49) in two steps: 

— For y < 2RM: On one hand, we have G(x) 4>{y) < C y. On the 
other hand, Hypotheses 1.7 and 1.3 show that b(y,x) > Cy 1-1 , so 
for y < 2RM (49) holds, as 7 < 2. 

— For 2RM < y < 2Rx: We have, again using (56), 

G(x)<j)(y) < Cye- a2X ~' < C y e'^ 1^ w) < Cyi~ x 

and wc conclude as in the previous case, by means of Hypotheses 1.7 
and 1.3. 

• Finally, considering £(y) = 1 and Hypothesis 1.7, we obtain C(z/) 1 = 
y^ 1 < C j/ 7 " 1 for any y > M because 7 > 0, and therefore (50) holds. 

□ 

Proof for growth-fragmentation. As in the self-similar fragmentation case, we 
only need to show that for growth-fragmentation model, K, M, R can be chosen 
appropriately so that equations (44), (45), (46), (49) and (50) hold for G, </>. In 
this way, as a direct application of Theorem 2.4, point 1 of Theorem 1.9 holds 
in that case. 

First, for the case b(x,y) = 2Bb/x with Bb > a constant, the bound 
(67) holds, and <fi(x) — for some constant > 0, as remarked at the 
beginning of section 4.2. In this simpler case, (44) is a consequence of (61), 
(45) is a consequence of (67), and (46) obviously holds with ((y) = 1, since 
4> is a constant. Similarly, (49) is obtained from (61), and (50) is true with 
£(y) = 1 and K = 2Bb- This allows us to apply Theorem 2.4 and prove point 1 
in Theorem 1.9 in this case. 

Let us consider now the case 7 > 0. Note that the requirements in eq. 
(64) hold, as due to (29) one may take /x = 1 in Hypothesis 1.4, and we have 
7 G (0, 2). Hence, all the bounds in Theorem 4.1 are valid here. 
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• The bound (44) is a immediate consequence of (61). With ((y) = y e , 
where < e < 1 and whatever M and R are, (46) is satisfied due to the 
fact that (j>(y) verifies (65): for Rz < y < 2Rz, 

7t4 < C(l + y 1 ^) < C(Co + C^) < C<j>(z). 

• Using (27), (65) and (66) we have for any x > M and for some constants 
denoted by C 

/•OO 

( / ) (y)G(y)dy<C (1 + y) dy 

Rx J Rx 

oo 

1+1 



<C y 1 e~ a ' dy = C e~ a ' < - Rx ^ +1 < C G(x), 

J Rx 



which holds by taking < a' < a 2 < (B m ) 2 /(B M (^ + 1)) and R > 1 such 
that 

a'R^ +1 > 



Bm 



7 + 1 
This proves (45). 

• As in the self-similar fragmentation case, we split the proof of (49) in two 
steps: 

— On the one hand for y < 2RM and x < y, using (65) and (62) we 
have G{x) <j>{y) < C y(l + y) and using Hypothesis 1.7 one sees (49) 
holds because < 7 < 2. 

— On the other hand for y > 2RM and x > ^ (which falls in the case 
max{2i?Af, 2Rx} = 2Rx) we have, again using (65) (66), 

G(x)cf>(y) < C (1 + y)e~ a ^ +1 <Cy>-\ 
and we conclude due to Hypothesis 1.7. 

• Finally, considering Q(y) = y t with < e < min{7, 1} and Hypothesis 1.7, 
we obtain ((y)y^ 1 = y £_1 < Cy 1 ~ 1 for any y > M, and therefore (50) 
holds. 

□ 



5.2 Spectral gap in L 2 space with polynomial weight 

Gathering the first point of Theorem 1.9 with some recent result obtained in [5] 
(see also [11] for the first results in that direction) we may enlarge the space in 
which the spectral gap holds and prove part 2 of Theorem 1.9 
We will make use of the following result 
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Theorem 5.1 ([5]). Consider two Hilbert spaces H and H. such that H cH and 
H is dense in T~L. Consider two unbounded closed operators with dense domain 
L on H , A on T~L such that A\jj = L. On H assume that 

1. There is G £ H such that LG — with \\G\\h = 1/ 

2. Defining </>(/) := (f,G) H G, the space H Q := {/ £ H; ^(/) = 0} is 
invariant under the action of L. 

3. L — a is dissipative on Hq for some a < 0, in the sense that 

VgeD(L)nH ((L-a)g,g) H <0, 
where D(L) denotes the domain of L in H. 
4- L generates a semigroup e tL on H; 
Assume furthermore on H that 

5. there exists a continuous linear form ^S> : TL —> R such that = ip; 
and A decomposes as A = A + B with 

6. A is a bounded operator from W to H ; 

7. B is a closed unbounded operator on H (with same domain as D(A) the 
domain of A) and satisfying the dissipation condition 

VgeD(A) ((B-a)g,g) n <0. 

Then, for any a £ (a, 0) there exists C a > 1 such that for any gi n £ H there 
holds: 

Vi > ||e tA g m - *(.g m ) G\\ H < C a e at \\g m - V(g in ) G\\ H . 

Proof of part 2 in Theorem 1.9. We split the proof into several steps. 

Step 1. Under the hypotheses of Theorem 1.9, take G, cf> solutions of (5), 
(6), respectively. We define H := L 2 (<j>G~ l dx) and H := L 2 (6dx) with 6 = 
<p(x) +x k , k > 1. Due to the bounds of G and (j> proved above, one can see that 

hch. 

We define 

A g := -a(x) d x g - (A + B(x)) g + C+g 
on T-L and L := A\h on H. We also define 

/>oo 

*(g) := (g,G) H = / g<t>dx (g £ U) 
Jo 

and ip := ^{h- From part 1 in Theorem 1.9 it is clear that L satisfies points 
1-4. Moreover, ^ is correctly defined and continuous on T-L as soon as k > 3, 
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so that H C L}{4>dx) thanks to Cauchy-Schwarz's inequality. To finish proving 
point 5, for given M, R > 0, we define x : = M l[o,i?,] i 



and 



-gx~ a(x)d x g - Xg 



C+g - 



= A(^) - gx, 



so that A = A + B and clearly A satisfies point 6. In order to conclude we have 
to establish that B satisfies point 7 for some well chosen k, M and R. Let us 
prove this separately for the cases a(x) = x and a{x) = 1. 

Step 2. The self-similar fragmentation equation. For a(x) — x, one has 4>(x) = x 
and we may easily compute, for m > 1, 



{Bg,g)m x 



with 



, dx) = Ti + T 2 + T 3 + T A + T s 



d x (x m x) 2 m + 1 
5 2 = — H~ 



T 2 :=J^(-2g)gx m = -2 J ' g 2 x m dx, 

£ B(x) gx m g< -B m J g 2 x m+ ^ 



g 2 x m dx, 



T 3 
T 5 



g 2 x m x , 



{£+g)x m g- 



Introducing the notation Q{x) — \g(y)\ y 1 1 dy, we compute, using (29), 



poo / f*C 

T 5 < J x m \g(x)\(p M B M J 

poo 

/ 2QQ' x m+1 -~< dx 
Jq 

POO 

■ g 2 d x (x m+1 -~<)dx. 

Jo 



\g{y)\y 7 1 dy)dx 



PmBm 
2 

PmBm 



Thanks to Cauchy-Schwarz inequality we have for any a > 1 

1— o poo 



Q 2 (x) < / y 2 ^ 2+a g 2 (y)dy / y~ a dy < 



a- 1 



y^-^g\y)dy. 



We then deduce 

PmB m m + 1 — 7 



2 a- 1 

/•GO 

< W y" +m g 2 (y)dy, 



2 ^- 2+a g 2 {y) / x^x™- 1 dxdy 
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with 

, v PmB m (m + 1) — 7 1 
v = i/(a, m) = fi(a, to), fi(a, to) := — — - x 



(to + 1) — 7 — (a — 1) a — 1 ' 

provided that to + 2 — 7 — a > 0. In particular, we notice that for m = 1 and 
a* = 2 - 7/2 e (1,2) we have f(a*,l) = 2P M B M /(2 - 7), so that for to = 1 



5 " 2-7 



/ x 1+1 g 2 dx =: C / x 1+1 g 2 dx. 
Jo Jo 



We also need to use the above calculation for to = k. We can hnd a and k such 
that 

k > 3, Ka<fc + 2-7, i/(a, A;) < 



To see this, take 



2 / ,! . / 1 1 \j 1 / 
a = 1 H , so that < 



B m a - I 2p m BM 

and then take k large enough so that 

fc + 1-7 < 2 
fc+l-7-a+l _ 

Putting together the preceding estimates we have proved 

/>oo 

(Bg,g) n < / x 5 2 (x) {-1 - X + (C - 1) x 7 } ^ 
Jo 

+ J™ ^g\ X ){ k —^- X -^ x ^dx. 

Recalling the definition of x, for an Y C > we can find R and M large enough 
so that 



poo 

(Bg,g)n < -C 9g 2 dx, 
Jo 



and that proves that assumption (7) in Theorem 5.1 is fulfilled with a = C = (3. 
The conclusion of Theorem 5.1 provides the conclusion in Theorem 1.9. 

Step 3. The growth-fragmentation equation. In this case we have a(x) = 1 and 
<j) is the solution to the dual eigenvalue problem (6). We first compute 

{B 9-,g)L 2 (<f>dx) = T123 +T4 + Is, 



31 



with 



p oo 

T123 ■= / {-d x g - Xg - B g} gt 
Jo 



1 f°° 

= -- J { + 5 </> + £+>} .g 2 <0, 

f°° 2 

Ta := - / g x<t>, 
Jo 

poo 

Tb ■= / {C + g)<t>g 
Jo 

< J™ B M Ci (1 + x) \g(x)\ (jf° | 5 (y)| y^ 1 dy) dz 

< C 2 / 3 2 a; 7 (1 + x) dx 

Jo 

poo 

< C 3 5 2 (0(x)+x 1+7 )dx, 



where as in the previous step we define Ci = (BmCi/2) (/x(a, 0) + fJ,(a, 1)) for 
some 1 < a < 2 — 7 (recall that here we have made the hypothesis 7 £ (0, 1)) and 
C3 comes from the fact that 4> is uniformly lower bounded by a positive constant. 
Following the computation of the previous step we easily get an estimate on 
(£> g,g)L 2 (x k dx), choosing k as before. Putting all together we obtain 

(Bg,g) n = [°° {B g) g {</>{x) + x k ) dx 



< / {C 3 ^(x)+x 1+ ^)-x<f>}g 2 dx 



{^x k -i-(\ + X )x k -^x k ^}g*dx. 



Again, for any C > we can find R and M large enough so that 

(Bg,g) n <-C 9 g 2 dx, 
Jo 

and we conclude as in the previous step. □ 



6 Appendix 

The following results are useful for dealing with weak conditions on the frag- 
mentation coefficient b(x,y). 
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Lemma 6.1. Let {fi}i<£i be a family of nonnegative finite measures on [0, 1], 
indexed in some set I , and take k > fixed. The following two statements are 
equivalent: 

Be,8e(0,l): [ ft < 6 f /< for all tel. (82) 

J[l-e,l] J [0,1] 

3 P G (0, 1) : / x k fi(x) dx<P I ft for all i G I. (83) 
J [o,i] J [0,1] 

Proof. First, assume (82) holds for some e, 6 G (0, 1). Observe that (82) is easily 
seen to be equivalent to 

/ fi >(!-<*)/ /< for all i e I. (84) 

Using this, 



x fi(x) dx = / x fi(x) dx + x fi(x)dx 

[0,1] J[0,l-e) y[l-e,l] 

< (l-e) fc / fi(x)dx+ f fi(x)dx 

J[0,l-e) J[l-e,l] 

f % (x) dx - (1 - (I - e) k ) [ fi(x)dx 

[0,1] J[0,l-e) 



< (l-(l-(l-e) k )(l-S)) / fi(x)dx 

= (S + (1 - 6)(1 - e) k ) f fi(x)dx, (85) 



[0,1] 

where (84) was used in the last step. This proves (83) with P := 5 + (1 — <5)(1 — 
e) k < 1. 

Now, let us prove (82) assuming (83) by contradiction. Pick e, 6 G (0,1), 
and take i G / such that (82) is contradicted for these e, (5. Then, 

x k fi(x)dx > / x k fi(x)dx 

[0,1] ./[l-e.l] 

> (1 - e) fc / fi(x) dx>(l- e) k 5 f fox) dx. (86) 

J[l-e,l] J [0,1] 

Hence, choosing S close to 1 and e close to gives an i G / such that (83) is 
contradicted. □ 

Lemma 6.2. Le£ 6e a family of nonnegative finite measures on [0, 1], 
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indexed in some set I . The following two statements are equivalent: 

V<5 > 3 e > : [ f l < S [ ft for alii el. (87) 

■'[1-6,1] -'[O.l] 

There exists a strictly decreasing function k i-4 Pfe , 
wit/i < p/j < 1 , lim pk = 

and / x k fi(x) dx < pk / /) /or a// i E I. 
J[o,i] J[o,i\ 

Proof. Let us first prove (88) assuming (87). Equation (85) holds here also, so 
x k U{x) dx<{5+(l- S)(l - e) k ) [ fi(x) dx 

[0,1] ^[0,1] 

Choosing 6 small enough, and then k large enough, one can take pt so that (88) 
holds. 

Now, let us prove the other implication by contradiction. Assume (88) does 
not hold, so there is some 8 > such that, for every e > 0, (88) fails at least for 
some j G /. With the same calculation as in (86), choosing e = 1 — (l/2) 1 / fe , we 
have that for every k > 1 there is some i € I such that 

x k f t {x)dx > ^ / fi(x)dx. 
[0,1] 1 J[o,i] 

This contradicts (88). □ 



Lemma 6.3. Consider a fragmentation coefficient b satisfying Hypothesis 1.1, 
1.2 and 1.5, and take < k < 1. For every S > there exists an e > such 
that 

f{\ — e)x px 

y k b(x,y)dy>{\-5) y k b(x,y)dy (x > 0). (89) 



/ 

Jo 



Proof. Equivalently, we need to prove that for every 6 > there exists e > 
such that 

y k b(x,y)dy<S y k b(x,y)dy (x > 0). (90) 

(l-e)x JO 

Using Hypothesis 1.5, take e > such that (28) holds with 8/k instead of S, 
where n is the one in Hypothesis 1.2. Then, 

y k b(x,y)dy < x k b{x,y) dy < x k - I b(x,y)dy 

(l-e)z J{l-e)x K Jo 

r-X f-X 

= x k 6B(x)=x k S -b(x,y)dy<5 y k b(x,y)dy. □ 
Jo x Jo 
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Corollary 6.4. Consider a fragmentation coefficient b satisfying Hypotheses 
1.1, 1.2 and 1.5. Then there exists a strictly decreasing function fc t-> p k for 
k > with linifc^+oo p k = 0, 

Pk >l for ke [0,1), Pl = l, 0<p k <lfork>l, (91) 

and such that 

( y k b{x,y)dy<p k x k B(x) (x > 0, k > 0). (92) 
Jo 

Also, for each < k < 1 f/iere exists p' k > 1 suc/i f/ia£ 

y fe 6(x, j/) <iy > p' k x k B(x) (x > 0, k e [0, 1)). (93) 



Proof. Apply Lemma 6.2 to the set of measures {f x } x >o given by 

f x (z):=b(x,xz) (z£[0, 1]), 

for which Hypothesis 1.5 gives precisely (87). Then, by a change of variables 
and using Hypothesis 1.2, (88) is exactly (92). 

For the second part, fix < k < 1. Applying Lemma 6.2 to the set of 
measures {z k b(x, xz)} x> q gives p' k so that (93) holds, as this set also satisfies 
(87) (by Lemma 6.3). □ 

Remark 6.5. One can omit Hypothesis 1.2 in the previous corollary and still get 
the result for fc > 1 by taking f x {z) := zb(x,xz) in the proof. 
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